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Abstract. The conditions for disk evaporation by elec- 
tron thermal conduction are examined, using a simpli- 
fied semi-analytical 1-D model. The model is based on 
t he m echanism proposed by Meyer & Meyer-Hofmeister 
( 1994 ) in which an advection dominated accretion flow 
evaporates the top layers from the underlying disk by ther- 
mal conduction. The evaporation rate is calculated as a 
function of the density of the advective flow, and an anal- 
ysis is made of the time scales and length scales of the 
dynamics of the advective flow. It is shown that evapora- 
tion can only completely destroy the disk if the conductive 
length scale is of the order of the radius. This implies that 
radial conduction is an essential factor in the evaporation 
process. The heat required for evaporation is in fact pro- 
duced at small radii and transported radially towards the 
evaporation region. 
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1. Introduction 

Ever since their theoretical rediscovery, advection dom- 
inated accretion flows (ADAFs, Abramowicz et al 



1995 



Narayan & Yi |!994|jl995[ Ichimaru |1977| ) have been widely 
regarded as the most likely source of Comptonized X- 
ray radiation observed from many X-ray binaries. On the 
basis of observational evidence (Lasota 1996| Narayan et 



al. 1997, Hameury et al. 1997 and references therein) it is 



believed that ADAFs form the inner part of the accretion 
disk system, while the outer part is formed by a standard 
Shakura-Sunyaev disk (SSD, Shakura & Sunyaev 1973| ). 
Despite the fact that such a bimodal disk geometry has 
alrea dy been propose d a lo ng time ago (Thorne & Price, 



1975; Shapiro, et al. 1976, henceforth SLE), no satisfac- 



tory theoretical explanation for these disk transitions has 
so far been found. 

It seems plausible that disk surface evaporation is re- 
sponsible for the transition. By some mechanism, the up- 
permost layers of the disk are heated up faster than radia- 
tive cooling can cool them down. The resulting hot 'vapor' 
forms a corona on top of the disk. Part of this vapor then 



accretes towards the central object, while another part 
moves outwards via a tran ssonic wi nd or breeze (Meyer 
& Meyer-Hofmeister 1994 , Liu et al. 
Turolla 1998 ). At a certain radius R eV ap 



1997| , Dullemond & 
the entire disk 
has been evaporated, and the 'corona' therefore becomes 
a true ADAF within this evaporation radius. 

It remains uncertain what drives the surface evapo- 
ration. Several mechanisms have been proposed. When 
no corona exists beforehand, one can show that the up- 
per layers of the SSD are unstable with respect to ther- 
mal perturbations (Shaviv 
Tschiipe 



Wehrse [1986| ; Hubeny [1989 
Kley |1993D . These layers will heat up m a run- 



away fashion, thereby causing the evaporation of the up 
per layers of the disk. Also, acoustic effects (Ickc 1976) 
and magneto-hydrodynamical fluctuations and instabili- 
ties may produce a corona, in a way similar to the for- 
mation of the solar corona (Galeev et al 



Pringle 1992) 



1979. Tout 



The most promising mechanism for disk evaporation is 
electr on thermal conduction (Meyer & Meyer-Hofmeister 
1994, henceforth MMH). A pre-existing hot 'corona' in- 



jects energy into an extremely thin layer on top of the 
underlying disk, turning this material into new hot coro- 
nal plasma. This mechanism is studied in this paper. 

In their paper, Meyer & Meyer-Hofmeister treat the 
evaporation problem in a 1-D vertical way, and obtain 
semi self-consistent models of disk evaporation. The goal 
of this paper is to re-examine this evaporation mechanism, 
but without the aim of finding fully consistent models. In 
fact, due to the complexity of the gas dynamics in the 
evaporation region, it is questionable whether it is possi- 
ble to model this process in a consistent way by use of 
1-D methods (Abramowicz et al. 1997). Instead, the goal 
of this paper is to study the structure of the boundary 
layer connecting the ADAF to the disk, and to study in a 
qualitative way the mechanism of evaporation, the condi- 
tions under which evaporation is possible, and the extent 
of validity of a 1-D vertical treatment of this evaporation 
process. To this end, a second order ordinary differential 
equation (ODE), containing only the most essential in- 
gredients for the description of the evaporation process, 
is derived. This equation is subsequently reduced to a di- 
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mensionless form containing two dimensionless constants. 
The resulting dimensionless equation is then solved and a 
simple relation between the two dimensionless constants 
is obtained by requiring a self-consistency condition to be 
fulfilled. Finally, by transforming the dimensionless solu- 
tions back to the real problem of disk evaporation, several 
important consequences are derived. 



2. Description of the model 

In order to explore the basic physics of evaporation by 
thermal conduction, we study the simplest possible sit- 
uation: a hot plasma flow on top of an accretion disk, 
with vertical thermal conduction, viscous heating, advec- 
tive cooling and radiative cooling. This simplified concep- 
tion of the problem should suffice for the goal of this paper. 
In principle one could think of this problem as the prob- 
lem of an ordinary ADAF and an SSD co-existing at each 
radius, and vertically glued together by a boundary layer. 
Such a model is consistent as long as the boundary layer is 
thin an d the evaporation is not too strong (Dullemond & 
Turolla 199§| ). Under these circumstances one can model 
the A DAF locally as a self-similar ADAF (Narayan & Yi 
1994). It is expected that the qualitative conclusions of 



the present model also remain valid somewhat beyond the 
breakdown of self-similarity. 

In the boundary layer the temperature drops steeply 
from the ADAF temperature down to the chromospheric 
disk temperature. The vertical scale associated with this 
gradient is extremely small, much smaller than the disk 
thickness. This justifies the use of a 1-D method for at 
least the lower parts of the boundary layer. The base of 
the boundary layer lies above the chromosphere of the 
disk, so that an optically thin treatment of the boundary 
layer is sufficient. The temperature gradient constitutes a 
heat flux pointing downwards towards the surface of the 
disk. At relatively high altitudes the ADAF produces an 
excess of energy which is transported downwards by the 
heat flux. In the lower parts of the boundary layer both 
radiative cooling and the upward gas motion absorb the 
flux. The upward gas motion, acts as a kind of vertical 
advective cooling (MMH). It is precisely large enough to 
cancel the downward flux, guaranteeing the flux at the 
base to be zero. This condition of zero flux at the base 
determines the evaporation rate. 

As the matter moves upwards from the disk surface, 
it has the tendency to shift towards larger radii. This is 
because, as the temperature increases, the pressure sup- 
port against gravity increases. A new gravitational force 
balance can only be reached by moving radially outwards 
a bit. But this phenomenon has no significant effect on 
the energy balance of the boundary layer, and we will ig- 
nore it. Also we will ignore the vertical friction (t z ^). We 
do take into account the decrease of f2 as matter moves 
upwards into the corona/ ADAF. 



The viscosity prescription is a delicate matter in 
boundary layers of the type studied here. Friction is as- 
sumed to arise from magneto-turbulence. The length scale 
associated with this turbulence is usually chosen ad-hoc, 
using the famous alpha-viscosity prescription: I = aH. 
For the disk height H one usually takes the rough esti- 
mate H — c s /VLki following from vertical pressure bal- 
ance. However, the length scales of magneto turbulence 
in the boundary layer are uncertain. If one would take 
I = ac s /Q.K then one finds that this length scale exceeds 
the boundary layer thickness, I 3> z , close to the lower 
boundary. This is inconsistent. In order to avoid having to 
discuss this highly uncertain issue, we rather ignore the ef- 
fects of vertical friction. In principle the radial friction t r ^ 
suffers from the same disease, but its contribution to the 
energy equation becomes small in the lower parts of the 
boundary layer, so this effect can be ignored here as well. 

3. The equations 

As a lower boundary to the calculational domain I take 
the height where the coronal temperature becomes equal 
to the chromospheric temperature. The height coordinate 
z is gauged to zero at this lower boundary, so that z = 
represents the height above the chromosphere. As an up- 
per boundary to the calculational domain I take a height 
z up obeying roughly z up < 0.3i? (R being the radius at 
which the corona is studied), so that geometric effects of 
the flow can safely be ignored. The corona extends well 
above this upper boundary, by virtue of the fact that it 
is assumed to be an advection dominated flow, for which 
the thickness H is roughly equal to the radius R. 

The model describes the temperature T, the density p 
and the vertical velocity v as a function of z. I presume 
that the system is quasi-stationary. The equations for T, p 
and v are the compressible Navier-Stokes equations in the 
coordinates R (radius), z (height) and <f) (azimuth), but in 
this calculation I choose the corona to be axisymmetric, 
and self-similar in the radial direction, thus reducing the 
problem to a one-dimensional problem in the coordinate 



The continuity equation is (pv)' — which integrates 



to 



pv = 
' 2 



(1) 



The radial motion and geometrical terms are neglected. It 
is assumed that the radial accretion predominantly takes 
place at altitude z > z up , allowing us to regard the evap- 
oration rate as a constant of motion in the domain of 
interest. The integration constant \& is the evaporation 
rate in units gem s . The factor 1/2 accounts for the 
fact that the disk has two sides. The pressure balance is 
(pc^Y = 0, where c s is the isothermal sound speed. This 
trivially integrates to 



P ci = u 



(2) 
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The ram pressure is neglected because the motion is as- 
sumed to be very subsonic. And as in the continuity equa- 
tion, the geometric terms are neglected here as well. 

The energy equation consists of five terms: viscous 
heating Q+, radial advective cooling Q a dv, optically thin 
radiative cooling <5_, electro-conductive heat flux J c and 
vertical advective heat flux J v . Viscous dissipation due to 
the vertical gradient of the rotational frequency, d£l(z)/ dz, 
is neglected. 

The dissipation due to the radial gradient in is given 

by 



0+ = pv (R 



L dR 



where the kinematic viscosity v is given by 

2 cl 

By using Eq. (||) the Q+ can be written as 

^ 3 2 
Q+ = 2 auJ 



(3) 



(4) 



(5) 



where the dimensionless rotational frequency lo is defined 
by = ujQk, with u> < 1. For an ADAF, the radial force 
balance equation, with omission of the radial kinetic term 
(which is very small), is 



n 2 R 2 JlOgW 



We take p oc i? _5//2 , so that this equation reduces to 

5 c 2 



lu z = 1- 



2Q 2 K R 2 



(6) 



(7) 



This can be substituted in Eq. (||), for use in the energy 
Eq. below. The radial advective cooling is given by 



de dp 1 



(8) 



where e = 0^/(7 — 1) is the thermal energy of the gas. The 
radial gas velocity is something that has to be estimated 
from the presumed radial structure of the corona. Take vr 
to be the usual expression 



n K R 



(9) 



If one wants to take into account the outwards shift result- 
ing from the requirement of pressure balance (discussed 
above), then one should add the velocity v^' defined as 

5 R v dc 2 



2 w 2 Vl 2 k R 2 dz 



(10) 



where v is again the vertical velocity. This extra velocity 
adds a contribution to the Q a dv, which is small enough not 
to influence the result significantly. For the sake of clarity 



it is ignored here, although it is easy to incorporate it. By 
using equations (0, ||), Eq. di) becomes 



7' 



1 



n 



n K R 2 



(ii) 



Radiative cooling is denoted with the symbol Q-. Sev- 
eral cooling mechanisms can can play a role, but the most 
important cooling mechanism in the boundary layer is 
Bremsstrahlung, Q_ ~ 5.0 x 10 20 p 2 - v / T. The tempera- 
ture T is related to c 2 by kT — pm p c 2 , where p is the 
molecular weight, and m p is the proton mass. The mean 
particle weights are p = 0.59, pi = 1.23 and p e = 1.14. 
The Bremsstrahlung cooling function is (using Eq. (||)) 



Q_ = K 1 n 2 c ; 3 

The constant K\ is defined as 

K 1 ~ 4.2 x 10 16 cm 6 erg _ V 
The vertical advective flux is J v , is 

Jv = pvh = \ — — -He 2 
2 7 — 1 



(12) 



(13) 



(14) 



where h is the enthalpy, and Eq. ([!]) has been used. The 
kinetic energy term has been neglected. Finally the con- 
ductive flux J c ~ -9.2 x 10~ 7 r]T 5 / 2 dT/dz (Braginskii 
^963] Spitzer ^962| ) can be written in the form, 



Jc 



-k c: 



dz 



(15) 



The conductivity coefficient Kq, for conduction along the 
magnetic field, is 



K ~ 2.8 x 10" 35 77 erg s 6 cm" 



(16) 



The Coulomb logarithm is roughly In A ~ 20. The coef- 
ficient r\ < 1 is put in as a fudge factor to parameterize 
the reduction in mean free path length as a result of pos- 
sible collective plasma modes and confinement by random 
magnetic fields. The microphysics of plasmas in these con- 
ditions is insufficiently understood to allow an estimate of 
77 from first principles, so we must retain it as the main 
unknown parameter of the model. 
The energy equation is now 



dJ c 

dz 



dJy 

dz 



— Q+ — Qadv — Q- 



(17) 



This is the basic equation of this paper. The density p and 
the velocity v have all been eliminated in favor of c s , and 
the equation has reduced to a single second order diffusion 
equation. 

The relation between the coronal accretion rate M c = 
—AirRHpvfi and the value II is an integral of the model 
over the entire vertical height H. But within the range of 
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validity of this model, a good estimate is given by assum- 
ing that the corona can be approximated as a homoge- 
neous flow of height H = \Jhj2 <Cs> / (where <Cs> is 
the average temperature in the corona). One obtains 



n : 



Mr n 



K 



27rx/T0 aRa 



(18) 



here the symbol a is defined as a ~<c^> /VIkR- It follows 
from the solution of Eq. ( f28|) below. For 7 = 1.5 and low 
enough M c this value is o ~ 0.59. 

4. Dimensionless form 

By defining the variable y as 



where £ is 



y 



£7/2 



3(7 

and a new coordinate x as 



2 3 
T) + 2 



21a 



Eq. ( fl7j ) becomes 



.fy +Dy -^ = 1 

ax z ax 



y 



2/7 



-3/7 



(19) 



(20) 



(21) 



(22) 



This is the dimensionless form of the main equation of 
this paper. From left to right one has the terms repre- 
senting the divergence of the conductive flux, the vertical 
advective cooling, the viscous heating, the radial advec- 
tive cooling and finally the radiative cooling. The symbols 
C and D are defined as 



C-- 



D-. 



2K 1 e /2 



n 



3a 



n 4 K R 3 



£3/2 $ 1 



7 - 1 V 2iaK vn n 2 R 3 / 2 



(23) 
(24) 



Using Eq. (pTSj) the constant C can be directly related to 
the coronal accretion rate M c by, 



C = 



3rr\/l0a 2 CT^i? 4 
The values of II and ^> are, 
3 a 



n 



* = 



2^3/2 c 



7-1 



63 K a 4 3 

IfVi^e 372 K 



CD 



(25) 

(26) 
(27) 



The solution for the corona at high altitude (x ^> 1) 
must be such that both fluxes J c and J v vanish. The equa- 
tion is then simply a balance between viscous heating and 
advective and radiative cooling, 



1 



y 



2/7 



Cy 



-3/7 



f) . 



(28) 



This equation is of rank 5 and cannot be solved exactly, 
but an very good parabolic approximation (to a few %) is 
given by 




28/11 



(29) 



For C < C cr it = 6vl5/125 this equation has two branches 
of solutions. The upper branch represents the ADAF 
branch y a daf-, while the lower branch represents the SLE 
branch y s i e (SLE 1976|) . For most of the solutions dis- 
cussed in this paper, the corona is found in the ADAF 
state, y(oo) = y a d a f- At C = C cr it = 6\/l5/125 the two 
branches meet. This value C cr it represents the critical ac- 
cretion rate for optically thin accretion flows (Abramowicz 
et al. |1995|), 



M c 



1871-76 
25 



« a p3 R 4 



(30) 



The M c andC relate to each other as M c / M cr i t = C/C c 



5. Solution without evaporation 

As a simple illustration let's first solve the equations with- 
out evaporation, so simply put D to zero. Eq. (|22|) becomes 



fy 

dx 2 



1 



y 



2/7 



Cy 



-3/7 



(31) 



x up the thermal conductive 



At the upper boundary x 
flux should vanishes: J c cx dy/dx = 0. At the lower bound 
ary x = the temperature should vanish, y = 0. The 
latter condition is of course unphysical, since the temper- 
ature should equal the chromospheric temperature rather 
than zero. But since the chromospheric temperature is pre- 
sumed to be very low compared to the coronal tempera- 
ture, this approximation is very good (within 1%). 

By taking the first integral of Eq. (|3l|), one finds the 
following expression for the dimensionless heat flux j c = 
dy/dx cx J c , 



dyV _ _ 2 
dx J 

+ 2yoo - 



14 



9/7 



-lev 



4/7 



where j/oo is the value of y(x — > 00), which is a solution 
to Eq. (p8[). This equation can be integrated numerically 
(using a relaxation scheme, for example), and the results 
are plotted in Figs, [j] and |2[ The result shows that at 
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x = x up the corona is in the ADAF state (y = y a daf)- 
Closer to the base of the corona the temperature drops and 
the conductive heat flux starts to grow. It has a maximum 
when y{x) = y s ie, at a certain x = x s i e . For x > x s i e there 
is excess heating and flux is being produced. For x < s i e 
there is excess cooling and the flux is being absorbed. At 
x = some flux remains unabsorbed. The value of dy/dx 
at x — can be found analytically from Eq. (||), 



dy_ 

dx 



2y 



14 9/7 
9 



(32) 



This non-zero flux at the base means that the corona is 
pumping energy into the chromosphere. But the chromo- 
spheric gas is not able to radiate it away quickly enough, 
so it must heat up. This inevitably leads to evaporation 
of the upper layers of the chromosphere. In order to make 
a more consistent model, the vertical motion should be 
taken into account from the start. 

6. Model with evaporation 

An upwards motion of the gas constitutes an additional 
vertical advective cooling (MMH). By allowing the con- 
stant D to be non-zero, this cooling takes effect. The value 
of D (being a constant over the entire domain) is deter- 
mined by adding an additional boundary condition to the 
system, 



dfj_ 

dx 



0. 



(33) 



x=0 



This yields both a solution for y(x) and a value for D. The 
evaporation rate is therefore found as an eigenvalue of the 
system, similar to the case of interstellar cloud evaporation 
s tudie d by Cowie & McKee ( |1977| ) and McKee & Cowie 
( |1977| ). The solution for the flux dy/dx is plotted in Fig. |. 
One can identify three regions in order of decreasing x, 

1. For x eq < x < x up the corona is nearly in local heat- 
ing/cooling balance. It is a solution of Eq. (|2"g|). 

2. For x m ax < x < x eq the heat flux rises from almost 
zero to a maximum at x ma x, as a result of an excess 
in local heat production. 

3. For < x < x max the heat flux decreases again to 
zero as a result of both radiative cooling and evapo- 
ration. For this region the solution is approximately 



a power-law, T 



oc y 



2/7 



Jr. 



-dy/dx oc x 2 / 5 . 



c 2 / 5 and the flux goes as 



The value of D, following from these models, depends 
on the value of C, or in other words on the coronal ac- 
cretion rate M c . This functional dependence of D on C is 
plotted in Fig. |]. The curve is closely fitted (to within a 
few %) by the expression 



(34) 



Note that, although the fit is very good, it is not exact. 
This expression then gives the evaporation rate D as a 
function of C. There are two branches. The upper branch 
has evaporation for small C (small M c ) and condensation 
for large C (large M c ). For this upper branch the corona 
(region |l|) is in the ADAF mode. The lower branch only 
has condensation, and the corona is in the SLE mode. 
Since SLE flows are known to be thermally unstable, this 
lower branch is not likely to represent any physical situ- 
ation. The dimensionless evaporation rate for the upper 
branch is then 



D 



3 

10 



2 


* /l- 12 > 






3 


3\/ 6VT5 



(35) 



The final formula for the evaporation rate is then (using 
Eqs. © and (HF 



£^5/2^ 
K 7= f{C/C cri t) 




(36) 



where f(C/C cr it) = (10/3) D(C) is the saturation func- 
tion, which is unity for C <C C' cr it, i.e. for M c <C M cr i t . 

There is a distinct value of C for which no evaporation 
nor condensation takes place, 



c - 21 c 



(37) 



which corresponds to M c — (21/25)M cr i t . For a more di- 
lute corona there will be evaporation, while for a denser 
corona there will be condensation. It is therefore tempt- 
ing to conclude that the corona will always tend towards 
saturation at M c = (21/25) M cr i t . But there is a caveat 
here, which will be discussed below. 

Before concluding this section it is interesting to com- 
pare the present model with the model of MMH. In the 
model of MMH part of the evaporated material leaves the 
system through a transsonic wind, while the other part ac- 
cretes radially onto the central object. By making a plau- 
sible assumption for the radial derivatives of the density 
and velocity, the evaporation problem is reduced to a 1-D 
vertical problem. The temperature in the wind is fixed by 
the conditions at the sonic point, while the pressure at the 
base, and thereby the evaporation rate, is determined by 
the balance between wind mass loss and conductive evap- 
oration rate. In this way a well-determined expression for 
the evaporation rate as a function of radius can be given. 

In the present model, a transsonic wind is not con- 
sidered, although it is not ruled out. The temperature at 
high z is fixed by the condition that viscous heating is 
balanced by radial advective cooling and radiative cool- 
ing. The evaporation rate is kept a function of the pres- 
sure, which is related to the radial accretion rate in the 
corona. Instead of making any assumption of the radial 
derivative of density, we have derived a relation between 
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the evaporation rate and the accretion rate in the corona. 
In the next section we will relate these to the the conduc- 
tive scale height and to a measure for the effectiveness of 
evaporation. 

7. Efficiency and conductive scale height 

The corona can only reach saturation (M c = (21/25)M cr it) 
when the evaporation is efficient enough. The evaporation 
rate has to compete with radial 'mass loss', i.e. the flow 
of coronal matter towards the central object. If a corona 
cannot evaporate the disk efficiently enough, the radial 
coronal flow will deplete the corona, until a low enough 
coronal density is reached for evaporation to compete with 
radial inflow. 

To investigate this one should compare the evaporation 
rate "J with the coronal accretion rate M c . Coronal mass 
conservation in a stationary situation is given by 



dM c 
~~d~R 



(38) 



Define an effectiveness index x as X = — d log M c /d log R. 
One has 



X = 2irR< 



Mr 



(39) 



The dimensionless number x gives the ratio between the 
evaporation rate and the coronal depletion rate due to 
radial accretion. One can also think of it as the ratio be- 
tween the time scales of vertical motion and radial motion. 
For x < 0.5 the evaporation rate is weak and for x > 0.5 
it is strong. It should be kept in mind that for \ > 0.5 
there e xist n o self-similar ADAF coronae (Dullemond & 
Turolla 1998 ), so that the evaporation models for \ > 0.5 
are not self-consistent with respect to angular-momentum 
conservation. By using Eqs.([39| [36], ^0|) one obtains 



3 / 7 c^KyKl 1-l jRg f{M c /M cnt ) 

R 



X 2 V 6V15 



7 



'M c /M c , 



X f(M c /M crU ) = 15.1^M f Jpi^d 

" lM c /M crit 



(40) 



where R g = 2GM/c 2 . In the last step j = 1.5. 

Closely related to the evaporation efficiency is the con- 
ductive scale height z c . This is the height below which the 
coronal energy content is significantly drained by thermal 
conduction. By using Eqs.([2l|, |30| ) one can find the dimen- 
sionless conductive scale height ( = z c /H, 



c 



= 7.2 




(41) 



The constant x c is the conductive scale height in the di- 
mensionless x-coordinate, introduced in section ^. As one 
can infer from the figures, the flux is roughly 20% of its 
maximum value at x ~ 4, so take the constant x c = 4. For 
(< 1 the thermal conduction only affects the lower layers 
of the corona, while for £ > 1 the conduction affects the 
entire corona. In fact, if Q > 1 one should expect radial 
thermal conduction to play an important role. 

If M c <C M cr n, the two dimensionless numbers x — X 
and f are roughly equal. For 7 = 1.5 and x c ~ 4 one finds 



C-1-9X 



(42) 



One sees that when the evaporation efficiency tends to 
become strong ( X > 0.5), the conductive scale height 
tends to exceed the size of the system £ > 1. This shows 
that for strong evaporation the conduction will dominate 
the entire corona from top to bottom, and radial ther- 
mal conduction will be an important effect. So, it is to 
be expected that for strong evaporation the problem be- 
comes essentially 2-D. A similar conclusion has already 
been drawn from angular momentum conservation con- 
siderations (Abramowicz et al. 1997 , Dullemond & Turolla 
|l998 ), but in that case a global radial 1-D solution with 
locally computed evaporation could still not be convinc- 
ingly excluded. The present conclusion is more dramatic, 
since it is unlikely that a very non-linear phenomenon like 
thermal conduction lends itself for dimensional splitting. 

Now consider the case when the corona becomes satu- 
rated, M c — > (21/25)M crit . Since M crit oc R^ 1 / 2 , this sat- 
uration can be only achieved when x = X f(M c / M cr i t ) = 
0.5, which implies x > 0.5 (since / obeys f(M c /M crit ) < 
1). This can be achieved for R < 10 3 (77/cv 2 ) Rg, but for 
most radii this means that x ^> 1 and therefore £ 3> 1. So 
one should conclude that coronal saturation is necessar- 
ily accompanied by strong radial conduction. Under those 
circumstances the present model breaks down, and fully 
2-D models should be used instead. 

8. Conclusion 

The model presented in this paper described qualitatively 
the mechanism of disk surface evaporation by thermal con- 
duction. It applies to accretion disc systems around black 
holes, neutron stars and white dwarfs. Although the model 
cannot be applied to the case of evaporation very close to 
a black hole (because of the relativistic velocities of the 
electrons, and the 2-temperature nature of the plasma) 
the qualitative picture sketched by this model is applica- 
ble to that case as well. 

In order for a disk to evaporate completely, the evapo- 
ration should be efficient. In a relatively short interval in 
radius, the complete disk should vanish, and the advec- 
tion dominated corona should be produced. The corona 
provides the energy for achieving this, and therefore the 
corona itself determines how strong the evaporation will 
be. In an equilibrium situation the evaporation rate should 
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balance the rate of radial inflow of coronal material. If 
most of the evaporation takes place in a relatively small 
range in radius, this means that the vertical velocity of 
the evaporating gas is in the order of the radial accretion 
velocity of the corona. The model described in this pa- 
per shows that in order to achieve this strong evaporation 
rate, the conductive scale height should be equal or larger 
than the radius. Consequently, radial thermal conduction 
will enter the problem. The heat required for evaporation 
does not come from the upper layers of the corona any- 
more, but is instead produced closer to the central object 
and radially transported to the evaporation region. There 
the flux will be directed down towards the disk surface 
and evaporate the disk. 

In order to model this complete disk evaporation with 
radial conduction, one should take the radial dimension 
of the problem into account. One could think of using the 
usual dimensional splitting procedure to solve the evapo- 
ration vertically (using a local 'heating' term to account 
for the divergence in radial heat flux), and the coronal 
dynamics and radial thermal conduction radially. Unfor- 
tunately this is procedure is highly questionable because 
of the very non-linear nature of the formula for conductive 
heat flux. It seems therefore that one should conclude that 
disk evaporation by thermal conduction is an essentially 
2-D process. 

These conclusions do not depend on the micro physical 
assumptions for evaporation, like the value of the conduc- 
tivity suppression factor rj. The conclusions here can be 
traced back to energy budget considerations. 

If electron thermal conduction is indeed responsible for 
the evaporation of a Shakura-Sunyaev disk, then there are 
also important consequences for the spectral modeling of 
ADAFs. If Comptonized soft-photons from the SSD con- 
stitute an important part of the spectrum, then this part 
emerges from the ADAF region that is strongly affected 
by radial thermal conduction. It is therefore question- 
able whether a consistent 2-phase (SSD+ADAF) spectral 
model can be built without 2-D radiative-hydrodynamic 
simulations. 
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Fig. 1. The temperature structure in dimensionless variables. Note that y 2 / 7 is proportional to the temperature. Both 
axes are logarithmic to clarify the structure. This solution is for D = and C = 0.15. 
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Fig. 2. The dimensionless conductive heat flux, plotted linearly. The solution is the same as for figure pi for D = 
and C = 0.15. 
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Fig. 3. The dimensionless conductive heat flux, plotted linearly. By requesting the flux at the base to vanish one finds 
D = 0.2055 for C = 0.15. 
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Fig. 4. The relation between the dimensionless evaporation rate D and the coronal accretion rate in the form of C. 



